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Abstract. The coordinate transformation between emission coordinates and inertial 
coordinates in Minkowski space-time is obtained for arbitrary configurations of the 
emitters. It appears that a positioning system always generates two different coordinate 
domains, namely, the front and the back emission coordinate domains. For both 
domains, the corresponding covariant expression of the transformation is explicitly 
given in terms of the emitter world-lines. This task requires the notion of orientation 
of an emitter configuration. The orientation is shown to be computable from the 
emission coordinates for the users of a 'central' region of the front emission coordinate 
domain. Other space-time regions associated with the emission coordinates are also 
outlined. 



PACS numbers: 04.20.-q, 45.20.Dd, 04.20. Cv, 95.10.Jk 
1. Introduction 

A relativistic positioning system consists of a set of four clocks A {A = 1,2,3,4) 
broadcasting their respective proper times by means of electromagnetic signals. 
Then, every event reached by the signals is naturally labelled by the four times {t^}: 
the emission coordinates of this event. Elsewhere [T], we have presented a brief report 
on relativistic positioning and related issues, providing a background with current 
references on the subject. 

A user of a positioning system that receives the four times {t^} knows his own 
coordinates in the emission system. Then, if he wants to know his position in another 
coordinate system, he must obtain the transformation between both coordinate systems. 
Thus, we must solve the following important problem in relativistic positioning. Suppose 
that the world-lines of the emitters 'jAi'T^) are known in a coordinate system {x"}: can 
the user obtain his coordinates in this system if he knows his emission coordinates {t"^}? 
Or, slightly more general, can the coordinate transformation = /t"(r^) be obtained? 
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The main goal of this paper is to solve this question in Minkowski space-time or, 
more explicitly, to obtain the coordinate transformation between emission and inertial 
coordinates for arbitrary world-lines of the emitters. 

In a two-dimensional approach to relativistic positioning systems this query is rather 
simple. Indeed, in this case the knowledge of the emitter's world-lines in a chosen null 
coordinate system gives the coordinate transformation between this null system and 
the emission coordinates at once [21 E]. This fact facilitates an analytical study of 
the positioning systems defined by inertial emitters in Minkowski plane [2], and those 
defined by stationary emitters in both Minkowski and Schwarzschild planes [3]. 

The general properties of the emission coordinates for the four-dimensional case 
have been analyzed in |3]. Nevertheless, in this generic case it is difficult to solve 
the above stated problem. The study of specific three-dimensional situations [5j has 
shed light on the complex geometry of the domains and codomains of the emission 
coordinates. On the other hand, the transformation between emission and inertial 
coordinates has recently been obtained for emitters following particular inertial motions 
in flat space-time, and also considering the immediate vicinity of a Fermi observer in 
the Schwarzschild geometry [B]. 

Here we present the solution to this problem for a generic configuration of the 
emitters in Minkowski space-time. We show that two emission domains exist which 
are called the front emission coordinate domain and the back emission coordinate 
domain, and we give the coordinate transformation for each one. The transformation 

= k°'{t^) between inertial and emission coordinates is given in a covariant 
way in terms of the world-lines 'jAi'r^) of the emitters. The compact covariant 
expression of our result is a powerful tool for use in subsequent applications: to 
obtain the coordinate transformation for specific configurations of the emitters or, under 
perturbation methods, modeling more realistic gravitational situations. 

The paper is structured as follows. In section [2} we pose the problem to be 
solved, defining the emission regions and the main relations governing the coordinate 
transformation: the null propagation conditions and the emission conditions. By 
choosing one of the emitters arbitrarily as reference emitter, in section [3| we show that 
the null propagation equations of a positioning system are equivalent to a rank-three 
linear system and a sole quadratic equation. Then, section |4] is devoted to obtaining the 
general solution of the linear system, and in section |5] we impose the remaining quadratic 
equation and we obtain a general form for both solutions. In section [6} we apply the 
emission conditions which guarantee that the solutions are physically interpretable as 
emission solutions. In section [7] we define the orientation of the positioning system 
with respect to an event, a concept which allows us to give the main result of this 
paper in compact form, namely, the explicit expression of the coordinate transformation 
between emission and inertial coordinates. The computational and/or observational 
determination of the orientation of the positioning system is analyzed in section [8j 
Section [9] deals with the analysis of our results and comments about ongoing work 
in progress and on further practical applications. Finally, an appendix is devoted to 
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presenting a technical proof of the resuhs given in section |4j 

A short communication on this work has been presented at the Spanish Relativity 
meeting ERE-2008 [7]. 

2. Statement of the problem. Emission regions and emission relations 
2.1. The emission region TZ of a positioning system 

Let us consider a positioning system in Minkowski space-time Ai^, and let 7a(t^), 
A = 1, ... ,4, he the world-lines of its four distinct emitters A (clocks) watching their 
proper time t^. These four times, broadcast by means of electromagnetic signals, will 
reach some region, say the emission region TZ, of events of Ai^. In the interior of TZ, 
the four times {r^} converging at every event define generically the so called emission 
coordinates [Zl El 13 El E] • 

It is clear that an event P of Ai^ belongs to the emission region TZ if and only if 
there exists a past directed null geodesic from P to every emitter 7yi(r'^) for some value 
of r^. For future discussions, we need to make the ingredients of this assertion more 
explicit. 

Let us denote by x the position vector of P with respect to the origin O of some 
inertial chart {x"}, x = OP, and by 7^ the position vector of the emitters with respect 
to this chart, 7^ = 07a(t"^). Then, in order for P to belong to TZ, the four vectors 

mA = X - 'JA , (1) 

which represent the trajectories followed by the electromagnetic signals in vacuum issued 
from the emitters A (see Fig. [T^), must verify the null propagation conditions L: 

L : (m^)2 = , VA . 

Furthermore, these four vectors have to be future-pointing or, in other words, must 
verify the emission conditions E : 

E : eu- niA < , VA , 

where u is any given, everywhere non vanishing, time-like vector field defining the arrow 
of time and 2e, e = ±1, is the metric signature. 

The events x where the emission condition E holds must be receivers. But it is 
worth noting that the null propagation conditions L also apply in the case of 'active' 
events x, able to send or to reflect null signals to one or more of the emitters. Such other 
location systems, as the physical realizations of coordinate systems are called here, have 
been dealt with elsewhere [TOl [TTl [12] but will not be considered here. 

Nevertheless, in obtaining our coordinate transformations, we will need to consider, 
besides the emission condition E, its 'causally dual' reception conditions R : 

R : eu ■ ruA > , WA , 

both condensed in the emission-reception conditions E-R : 

E-R : emA-ruB <0 , 'iA,B , 
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as it is easy to argue. 

Now, the above assertion about P and TZ may be stated as follows: an event P 
belongs to the emission region IZ of a positioning system of emitters 7yi(T'^) if, and only 
if, its position vector x in some inertial chart verifies the null and emission conditions 
L and E, respectively, for some values of the proper times of the emitters. 

2.2. The characteristic emission function Q of a positioning system and the emission 
coordinate region C of the space-time 

It then turns out that a positioning system may be considered as a device G that 
physically associates, to every event of the emission region IZ, a set of four times {t^}. 

From a formal point of view, this device is nothing but an application : TZ — > T, 
henceforth called the characteristic emission function of the positioning system, that, 
with every event x of the emission region 7Z of the space-time, associates four times 
of the grid T of the r's, : x i — > (r^) = G(x). 

4 

The grid T is nothing but the Cartesian product T = x {r} ^ M of the spaces 
(real lines) of the variables r (for details on the concept and role of the grid, see for 
example HEUS]). 

In this grid T, the image 0(7^) of the emission region TZ by the characteristic 
emission function G, is called the emission co-region of the positioning system, and is 
denoted by ^ = G(7^). The points of this region of the grid T are the quadruplets 
of times that can really be received in the space-time, so that they are related to the 
space-time events at which they are measured. In this sense, this emission co-region 
is the sole region of the grid T which possesses a physical meaning (the other quadruplets 
of T are a convenient mathematical completion of the emission co-region but with no 
relation to space-time events). 

Obviously, if the world-lines of the four emitters are sufficiently smooth and 
broadcast their proper time continuously, the emission co-region ^7Z C T is connected 
and, because of the regularity of the light cones in Minkowski space-time, TZ is connected 
too. But this property does not guarantee that the emission function G, = G^(x), 
is invertible in TZ. 

For G to be invertible, the gradients dr"^, normal to the hypersurfaces = 
constant, have to be well defined and linearly independent. But, because in Minkowski 
space-time our light cones are everywhere differentiable up to on their vertices, i.e. on 
the emitter world-lines 7a, the gradients dr^ are well defined everywhere in TZ up to on 
T^n (UA7A(r^))- Then, on the region TZ — (Uyi7yi(T^)) , because the dr"^ are metrically 
collinear to the above null vectors m^, the linear independence of the dr"^ may be 
expressed by the coordinate condition C: 

C : nil A A A ^ . 

This condition C is equivalent to say that jQ{x) ^ 0, where je{x) is the determinant of 
the Jacobian matrix Jq{x) of G, so that the locus where the dr^ are linearly dependent 
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is the hypersurface ^7 in 7^ of equation J' = {^0(3;) = 0}. We shall call the regions 
V = J U (Ua7a(t'^)) and C = 1Z — V the emission degenerate region, and the emission 
coordinate region respectively, TZ = C WD. 

It is to be noted that, the condition C being satisfied on its events, C is an open set. 
For the same reason, is invertible in C, so that all its events may be locally labeled 
by the coordinates {t"^}. Nevertheless, this does not means necessarily that C be a 
coordinate domain of a local chart (C, G), because the condition C assures only the local 
invertibility of G. In fact, we will see that C is not a coordinate domain, but the union 
of two coordinate domains; this is why we have called C the emission coordinate region. 
The region in the grid T where the characteristic emission function G, r"^ = G'^(x), is 
locally invertible is the region ®C = G(C), and will be called the emission coordinate 
co-region of the grid T. 

2. 3. The main relations of a positioning system 

When the world-lines of the emitters are known, the main relations of a positioning 
system are the set {L, E} of the null propagation conditions L and the emission conditions 
E. 

In terms of the position vectors 7^1 of these world-lines '~^a{j^)) they give rise, by 
([1]), to the null propagation equations 

L: (x-7a)-(x-7a) = , MA. (2) 

and to the emission inequalities 

E: eM-(x-7A)<0, VA , (3) 

where u is any given, everywhere non vanishing, future-pointing time-like vector. 

To invert the function G in the emission coordinate co-region ®C is to solve the 
mam relations (g, Q in X, X = ^(r'^), for values of the r^'s verifying the coordinate 
condition C or, by ([T]), the coordinate inequality: 

C : (x - 71) A (x - 72) A (x - 73) A (x - 74) ^ . (4) 

These solutions x = h{t^) may alternatively be read as the position vectors of the 
events whose past light cone cuts the emitter world-lines 7a(t'^) at their times t^. The 
components x" = k°^{t^) of x then define the coordinate transformation between the 
emission coordinates {t^} and the inertial ones {x""}. 

The main object of this paper is to obtain this coordinate transformation. 

3. The null propagation equations L 

When {t^} are the emission coordinates of the event x, the emitters are at the events 
{1a{.'t^)}- These four events define the internal configuration of the emitters for the 
event x. 

One can try to solve the null propagation equations ([2]) in x, x = ^(r^) directly, 
but it is better to first carefully separate the ingredients intrinsically related to the 
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Figure 1. (a) The light-like vectors to^ = x — '^a represent the trajectories followed 
by the electromagnetic signals between the emitters and an event P in the emission 
region, x is the position vector of this event with respect to some inertial chart, and 
7^1, A = 1,2,3,4, are the position vectors of the emitters with respect to this chart, 
(b) If we choose the emitter 4 as origin (reference emitter), the relative positions of 
the others (referred emitters) are = ja ~ 74i ^ = Ij 2, 3, and the position vector of 
the event P is y = m^. 



configuration of tlie four emitters A, which are independent of the inertial chart {x"}, 
from the ingredients related to this chart, which is independent of the configuration of 
the emitters. 

For this purpose, we shall arbitrarily choose one of the emitters, say A = 4, as 
reference emitter and we shall search the solution x under the form: 

a; = 74 + 1/ , (5) 

where y is the solution to the null propagation equations ^ when the origin is chosen 
at the position 74 of the emitter 4 (see Fig. Thus, y is submitted to 

L4: {y-eA)-iy-eA) = , VA , (6) 

where 

ca = 7a - 74 • (7) 
Now, splitting these relations for A = 4 and A = a = 1, 2, 3, we have 

64 = , eA = ea (8) 
and, for the null propagation equations ([6]), 

y' = o, A = 4, 

y2_2ea-y + iea)' = 0, A = a. ^ ^ 

Noting that are the relative position vectors of the referred emitter A = a with respect 
to the reference emitter A = 4, and that the half of their squares, represent the world 
functions [13] of the referred and reference emitters, we have: 

Lemma 1 The null propagation equations ^ are equivalent to the main linear system 
Ca-y = fla , a = 1,2,3 , (10) 
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and the main quadratic equation 

y' = 0, (11) 
where the scalars VLa, world functions of the endpoints of the vectors Ca, are given by 

= \{eaf . (12) 

In the study of the null propagation equations we are supposing that the emitted 
times {t^} received at x constitute the emission coordinates of x so that the coordinate 
inequality Q fulfils. From the definitions ([T]) this condition writes 

C4: X- 1/7^0, (13) 

where x is the configuration vector defined by 

X = *(ei A 62 A 63) , (14) 

with * being the Hodge dual operator. Note that, under the form C4 of C, x is necessarily 
non vanishing: 

Xt^O, (15) 

which means that the configuration {7a(t^)} determines a hyperplane. This 
configuration hyperplane contains the four events {ia^'t^)} and is orthogonal to the 
configuration vector x- The sign of x^ provides the causal character (space-like, light- 
like or time-like) of this hyperplane. Thus, in the emission coordinate region C one can 
distinguish three disjoint regions: 

the space-like configuration region Cs = {x&C\tx^< 0}, 
the null configuration region Ci = {x E C \ = 0}? 

the time-like configuration region Ct = {x&C\tx^> 0}. 

with C = U C£ U Ct- 

The ingredients Via and x of the internal configuration of the emitters make use of 
the relative position vectors Ca and their double exterior product ci A 62 A 63. Let us 
complete them by introducing their simple exterior products in the form 

= *(ea+i Aea+2) , (16) 

where the operations on the indices obviously have to be understood modulo 3. 

In looking for the solutions to the null propagation equations (|6]), it is convenient 
to first start solving the main linear system ((To|) in y and then constraint it to verify the 



main quadratic equation (11). Then, the solution(s) y = t{T^) so obtained, incremented 
by the position 74 of the reference emitter A = A, will give by ([s]) the solution(s) to 
the null propagation equations (|2]), x = k{t^) = 74 + ^(t^), i.e. the wanted coordinate 
transformation between emission coordinates and inertial ones. The next section is 
devoted to solving the main linear system (10). 
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4. The main linear system 



The regularity condition (15) for the configuration vector x means that the rank of the 
main linear system ( 10 ) is exactly 3. As a consequence, its general solution y depends on 
a sole parameter, say A, so that denoting by a particular solution, one has y = y* + Xx 
because, by (14), x is a vector such that Ca ■ x = 0. We can thus state: 

Lemma 2 In the emission coordinate region C, the general solution to the main linear 
system ( 10 ) is of the form 

+ Ax , 



y 



(17) 



where the parameter A takes arbitrary values, x is the configuration vector ( 14 ) and y^ 
is a particular solution. 

A tentative study of explicit expressions of the particular solution y^. shows that, 
for every emission coordinate domain, there is not a sole analytical function, depending 
exclusively of the elements of the configuration {7a(t'^)} of the emitters, that is valid 
throughout the domain. This is why it is necessary to introduce an external element to 
obtain a sole analytical expression. We shall see that it is sufficient for this external 
element to be a vector field ^ transversal to the configuration, i.e. such that 



and otherwise arbitrary. Denoting by i{ 



the interior product, it may then be shown 



see 



Appendix A for the proof): 



Proposition 1 In all the emission coordinate region C, the general solution to the main 
linear system (10), 

Ca - y = ^a , a = 1,2,3 , 

is of the form 

y = y, + Xx, (19) 
where the parameter A takes arbitrary values, x is the configuration vector 
X = *(ei A 62 A 63) , 

and y^ is the particular solution orthogonal to a chosen transversal vector ■C; ^ " X 7^ 0, 
given by 



1 



i ■ X 



(20) 



where the bivector H is, like the vector x, ^ function of the configuration of the emitters, 

E'' = *{ea+i A ea+2) . (21) 



H = VtaE" 
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The direct and simple discussion that has enabled us to reach proposition [T] is the 
result of a careful search in order to find, in every coordinate domain, a sole, general 
and covariant expression for the coordinate transformation between emission coordinates 
{r"^} and inertial ones {x"}. 

The last two requirements of generality and covariance enable us to find once and 
for all and for any inertial system (and in compact form due to the intrinsic vector 
formalism) the coordinate transformation in question. But they do not suffice to lead 
by themselves to a sole expression valid for all the configurations, i.e. valid at all the 
events of every emission coordinate domain. For example, if we were to work directly 
on the main linear system in the cases of a configuration generating a nonnull vector x, 
we would find the solution 

y = \z{x)H + Xx, (22) 
X 



as now directly follows from (20) and (19) for the choice ^ = X- But, meanwhile in fact 
there is neither physical nor mathematical reason for the solution to lose analyticity, 
this expression (22) becomes manifestly undefined on the null configurations, = 0- 
This is why it is necessary to introduce an external element to find the sole regularized 
analytical expression above mentioned. 

This external element, the transversal vector field ^, independent of the confi- 
guration of the emitters, is otherwise arbitrary. But, meanwhile due to lemma [2] we 
know that the main linear system depends on a sole parameter A, the expressions (19), 
(20) of y apparently seem to indicate that it depends on the (4 + 1 =) 5 arbitrary 
parameters A}. Thus, what is the precise role played by this transversal field ^? 

From the definition (21) of the bivector H it is clear that its dual *H is orthogonal 
to X) "^(x) * -f^ = 0, so that, one has H = x ^ cl for some vector a. With this expression 
of H it is easy to prove from (19) and (20) that the variation of the solution y with 
respect to the parameters ^ verifies 

|| A X = , (23) 

which shows that the changes in ^ may be absorbed by A, i.e. that the parameters 
are not essential in the Eisenhart sense [Hj. Nevertheless their elimination gives rise 
to different expressions for the causally different regions that one can find in the whole 
emission coordinate region. 



5. The main quadratic equation 

Obtaining the transformation = K"(r^) between emission coordinates {t^} and 
inertial ones {x"} amounts to determine the intersection of the four future light cones 
emitted by the emitters at the configuration {7a ('T"'^)} or, in the dual interpretation, 
to determine the vertex of the past light cone that contains the emitters in such a 
configuration. 
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According to proposition [T| the solutions y of the main hnear system ( 10 ) depend 
hnearly of the parameter A, and consequently, describe a straight line in Ai^. The 
vertex that we are searching for is consequently on this straight line and, according to 
our choice of reference emitter, is the point on it related to the emitter A = Ahj the null 
vector 7714 = X — 74 = This is nothing more than what the main quadratic equation 
(11) expresses: y"^ = 0. This cone, whose vertex is the emitter A = 4, cuts the straight 
line ( 19 ) at a value of A given implicitly by: 

x'X' + 2{x-y.)X + y',=0 , (24) 

and it is this value of A which, put in the expression (19) of proposition [l| enables us 
to obtain the coordinate transformation y = l{t^) that we are searching for. 

But, before obtaining it, and in order to be sure of its real character, it is convenient 
to observe that the sole assumption we have made so far to solve the null propagation 
equations is that of the regularity of the configuration data {7^1 (t"^)}, i.e. of the non- 
vanishing of the configuration vector X 7^ 0. This assumption is only a strict part of 
the coordinate condition C, as its form C4, X • y 7^ 0, shows. 

The study of degenerate configuration data x = and of the degenerate coordinate 
condition x ■ ?/ = is interesting to gain more in-depth knowledge of the physical 
conditions leading to their presence, in order to avoid, control or predict them. 
But we shall restrict ourselves here, as already stated, to obtaining the coordinate 
transformation between emission coordinates {r"^} and inertial ones {x"}. For this 
purpose, from now on, we shall work under the coordinate condition C. 

In fact, it is this condition that guarantees the real character of the solutions of the 
null propagation system. To demonstrate this, we shall separately analyze the cases of 
null and non null configurations. 

Null configurations. In the region where the configuration vector is null, = 0) 
expression (19) for the solution y, y = y* + Ax, leads to x ■ 1/ = X ' Z/*, so that the last 
expression (13) for the coordinate condition C ensures the well defined character of A in 
(0, 

A = — -^2— , (25) 
2(x-y*) ^ ^ 

and we have 

Proposition 2 Under the coordinate condition C, for null configurations, x^ = 0, the 
null propagation equations (|6| admit a real and single solution y given by: 



where y^, and x o-^e respectively given by (20) and (14). 

Non-null configurations. In the regions where the configuration vector is non-null, 
XV 0, (|23) gives: 

A± = p (-X ■ 2/* ± v^) , A = (x ■ y.f - x^yl • (27) 
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But, for these non-null configurations, we can choose ^ = x the expression (20) 



of the particular solution to the main linear system (10). Let us denote by c the 
corresponding y^: : 

c=-^i{x)H. (28) 
X 

Because H is antisymmetric, it is obvious that c is orthogonal to x, so that, for the 



corresponding value Ac of A, (27) takes the simple form: 

/a: 



A 



±- 



X 



A, 



2 2 

-X c 



(29) 



Now, expression (19) for the solution ?/,?/ = c + Ac leads, by product by x, to ?/ ■ x 



AcX^. Taking its square and substituting in it the value (29) of Ac, one obtains: 

{yxf = ^c , 



(30) 



so that the form (13) of the coordinate condition C ensures the strict positiveness of the 



discriminant Ac- Then, for y^ = c and the values (29) of Ac, equations (19) lead, after 
a little analysis, to the following result: 

Proposition 3 Under the coordinate condition C, for non-null configurations x^ 7^ 0, 
the null propagation equations (|6| admit two real solutions y± given by 

X2 

y±=c±\c\u, u=j-^x, (31) 



where c and x are, respectively, given by (28) and (14). 



In the above proposition and in what follows we denote \v\ the modulus of a vector v, 

c, the discriminant A given by 



|t>| = y |f 

We have seen that, for the particular solution y^, 



(27) takes the value Ac given by (29). To what extent the discriminant A, corresponding 
to other particular solutions y^, ^ c, differs from the value Ac? 



The question is pertinent because we want to unify expressions (26) and (31) so 



as to have a single expression for every coordinate domain, irrespective of the causal 
orientation of the configuration vector x- To answer this, note that c is the particular 



solution to the linear system orthogonal to x and that, from (19), any y^: differs from 



c by a term of the form Ax- More precisely, c 
calculation leads to: 



y* 



'-X- Then, a straightforward 



Lemma 3 Under the coordinate condition C, for non-null configurations x^ 7^ 0, the 



discriminant A and each one of the solutions y± of the quadratic equation (24), given in 



(27), are independent of the particular solution chosen and, therefore, of the subsidiary 



vector ^ too. The discriminant is positive and its invariant value is 
A = Ac = -xV>0 , 



(32) 



where c is given by (28) 
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Unified expression for any configuration. As already stated, we want to find, 
for every coordinate domain, a sole expression valid for the whole domain, i.e. valid 
for all the configurations {7a (t^)} that could correspond to the events of the domain. 
Nevertheless, in principle, by continuously changing of events in the coordinate domain, 
we can make these configurations to have a vanishing, positive or negative value of x^- 
When this is the case, none of the two standard expressions (27) of A for 7^ comes 
down to the sole solution (25) for = ^ind, consequently, none of the expressions 
(31) reduces to expression (26), both of them becoming degenerate. To correct this 
degeneration, it is sufficient to multiply numerator and denominator of the second 
member of (27) by the conjugate of the numerator. Once this is done, taking into 
account proposition [T] one has the following result. 

Theorem 1 In all the emission coordinate region C, the solutions to the null 
propagation equations ^ are real and admit the expression: 



y± = y* + A±x , A± 



ix ■y*)±y/A 



^ = {x- y*f - x^y 



2 

* 1 



(33) 



where ?/* and x ore respectively given by (20) and (14). For null configurations x^ = 0, 
the sole solution is given by one of these expressions, the other becoming degenerate. 



6. The emission conditions E 

The expressions (33) of theorem [l] give the Minkowski events y that can be related to the 
configuration events {7a (t^^)} by means of emission or reception of null signals. But we 
want y to be the events reached by null signals emitted from the configuration events. 
Thus we have to impose on y the emission conditions E. This will select a special class 
of configurations, which will be called emission configurations. 

In a first step we shall impose the intermediate emission-reception condition E-R, 
guaranteeing that either all the signals at y have been emitted by the configuration, or 
all the signals at y will be received by the configuration; such a configuration will be 
called an emission-reception configuration. And in a second step we shall choose with E 
the first of these two cases. 

As we have seen in section [2| the emission-reception condition E-R states that all 
the null vectors joining 7^ and y, rriA = y — 'jA, must have the same (past or future) 
orientation, i.e. em a ■ ttib < ,^ A,B . In terms of the relative position vectors Ca of 
the configuration {7a (t^^)} given by ^ and (|8|, one has: 

"^4 - "^a = 7a - 74 = ^^^^ 

rria - rUb = - la = eb - Ca , 
so that, because {rriA — rnsY — — 2m^ ■ m^, we have: 

Proposition 4 For a regular configuration {lAij^)} to be an emission-reception 
configuration it is necessary and sufficient that all their relative positions be space-like. 
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In other words: the null directions tua verify the emission-reception condition E-R if, 
and only iff, 

(E-R)^: e(e,)2>0 , e (e, - efe)^ > . (35) 

Let us note that, in particular, this proposition tells us that, in the grid T ~ M"^ 
of the r's, the domains of emission coordinates {t^} are in the interior of the region 
determined by (35), so that the points of the grid T in the complementary region 
certainly have no physical meaning. 

Equations (33) of theorem [l] show us that, under the coordinate condition C, two 
real and definite solutions to the null propagation system may correspond to every 
configuration {7^1 (t^^)}- Now, under the additional emission- reception condition E-R, 
one of these definite solutions may be that of emission and the other of reception. This 
last one finally has to be detected and discarded by means of the emission condition E, 
which, under the E-R one, reduces to: 

E4 : e?/-M<0, (36) 

where u is any given, everywhere non vanishing, future-pointing time-like vector. 

For non-null configurations, 7^ 0, the general expression (33) gives the two 
admissible solutions under the coordinate condition C. These solutions also admit the 
non-unified expression given in proposition [3j From it we obtain: 

y+-y-=2^ (37) 



Then, (32) implies that if ex > {respectively, ex < 0) then ec < {respectively, 
ec^ > ) and, as a consequence of (37), ?/+ is future-pointing iff ?/_ is future-pointing 
{respectively, past-pointing). Thus, we have: 

Lemma 4 If x is a time-like vector, ex^ < 0, then only one of the solutions (33) 
corresponds to an emission configuration. 

If X is a space-like vector, ex^ > 0, then the two solutions of (33) correspond to 
either two emission configurations or two reception configurations. 

When X is a time-like vector, can we detect the solution which corresponds to an 
emission configuration? The answer is affirmative. Indeed, as a consequence of lemma [3} 
each one of the solutions y± is independent of C,- Thus, taking for them the non-unified 
expression given in proposition^ we obtain (when ex^ < 0): 



e(ex) ■y± = x-y± = ±|c| Ixl , 

and, consequently, we can state: 

Proposition 5 For a space-like emission configuration, ex^ < 0, the sole emission 
solution is y+ (respectively, y-) if ex is past-pointing (respectively, future-pointing). 

When X is a space-like vector, we can know a priori whether the two solutions y± 
correspond to emission configurations. Indeed, in this case c is a time-like vector. Then, 
making use again of the non-unified expression for y± given in proposition |3| we obtain: 

e c ■ y± = ec^ < , 
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and, consequently, we can state: 

Proposition 6 For a time-like emission configuration, ex^ > the solutions y± 
correspond to emission (respectively, reception) configurations if c is future- pointing 
(respectively, past-pointing). 

Finally, for a null configuration, = 0, proposition |2] gives the sole admissible 
solution under the coordinate condition C, which of course also reduces to one of the 
solutions of the general expression (33). In this case A = {y^ ■ x)^ and, consequently, 
the non degenerate solution in (33) is y+ {respectively, yJ) if e (ex) ■ y* = X ' U* > ^ 
{respectively, < 0). Thus, we can state: 

Proposition 7 For a null emission configuration, = 0, the non degenerate solution 
is y+ (respectively, y^) if tx is past-pointing (respectively, future-pointing). 

It is worth remarking the attachment between the results stated in propositions |5] and 
[7| If we continuously change to a null configuration coming from a space-like one, 
the configuration vector x keeps its future or past orientation. If x is past-pointing 
{respectively, future-pointing) we must take the (continuous) expression yj^ {respectively, 
y-) in both the null and space-like regions. 



7. The coordinate transformation from emission to inertial coordinates 

7.1. Front and hack emission coordinate domains 

Theorem [T] shows that the emission coordinate region C is mapped with two local charts 
and gives analytical expressions for the transformation between emission and inertial 
coordinates for the two corresponding coordinate domains. 

These expressions of the coordinate transformation and the study of the emission 
conditions in section [6] show that the causal character of the configuration of the emitters 
differs for the two emission coordinate domains. Space-like and null configurations 
only admit one of the solutions given in theorem [T] as emission solution. And, of 
course, the coordinate domain of this solution contains by continuity a time-like 
configuration region. The coordinate domain of the other solution only contains a 
time-like configuration region. 

We see that the two coordinate domains differ enough in the causal character of 
their emitter configurations. We shall call front emission coordinate domain the 
coordinate domain that contains events with the three possible causal configurations. 
More precisely, contains all the space-like configuration region Cg, all the null 
configuration region and a part of the time-like configuration region Ct- We 
shall call hack emission coordinate domain the other coordinate domain, that only 
contains events with time-like configurations]!] In fact, = Ct — Cf . 

I Let us remember that it is usual to call coordinate domain the open set U of any local chart {U, (p) 
of the atlas defining a differentiable manifold. This appellation requires attention because a coordinate 
domain is not necessarily a domain, but simply a (not necessarily connected) topological open set. In 
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Moreover, the two time-like regions have the same codomain in the grid, namely, 
Q{C[) = Q{C^). This relation is very important. It says that, whatever be the four 
values {t^} received by a user in the coordinate region C[ (resp. C^), another user 
in the coordinate region (resp. C[) may receive the same values {t"^}. In other 
words: a user that only receives four times {t^} defining a time-like configuration of 
the emitters is unable to detect in what part of the region Ct, in C[ or in C^, he is. 

7.2. Orientation of a positioning system 

To be able to detect in what of these regions a user is, a notion of orientation of a 
positioning system is necessary. We give the following one: 

Definition 1 The orientation of a positioning system with respect to an event of its 
emission coordinate region C is the sign e of the coordinate condition scalar: 

e = sgn * (rrii Am2 Am^ A m^) . (38) 



We have seen in section 2.2 that, in the emission region TZ, the hypersurface JT" 
separates the emission coordinate region C in two open sets. We can now identify them 
with the above front and back coordinate domains and C^, respectively. Because JT" 
is the hypersurface where the coordinate condition C is not verified (vanishing Jacobian 
je{x) = 0), the non vanishing member of its inequality has a constant sign in every 
coordinate domain and , so that we have the simple but important result: 

Proposition 8 The orientation e of a positioning system is constant in every one of 
the coordinate domains and . 



The same way that leads to the form (13) of the coordinate condition C shows that 



e = sgn {y ■ x)- Because y is necessarily future-pointing for a emitted signal, if x is 
time-like or null we can be sure that the sign of {y ■ x) is the same as that of the sign 
of {u ■ x) for any everywhere non vanishing future-pointing time-like vector u. This last 
sign is plus or minus according to the past- or future-pointing character of ex- Taking 
into account propositions [5] and [7| one has the following: 

Proposition 9 In the regions Cg and Ce, the orientation e of a positioning system is 
given by 

e = sgn {u ■ x) (39) 
for any future-pointing time-like vector u. 

fact, meanwhile the front emission coordinate domain is genericaUy connected, the back emission 
coordinate domain is generically the disjoint union of four connected components [5]. 
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7. 3. Explicit expression of the coordinate transformation from emission to inertial 
coordinates 

We shall comment below in section |8] that for users in the region Ct the determination 
of the orientation needs additional information. But for the moment, Theorem [T] and 
the above propositions lead us to the following result. 

Theorem 2 Let •ja be the position vectors of the world-line equations 7yi(T^) of the 
four emitters of a positioning system with respect to an inertial coordinate system 
and {t^} their emission coordinates. In all the emission coordinate region C = C^UC^ , 
the coordinate transformation x = /t(r^) is given by: 

2 

x = 74 + (40) 



where ?/* is the quantity given by (20), x ^■^ configuration vector (14) and e is the 
orientation (38) of the positioning system with respect to the event that receives the data 

This is the main result reached in this paper. It rest to analyze now to what extent 



this expression (40) may be determined by a user from the data received by him. 



8. The central region of a positioning system. Computational and 
observational orientation 

The world-lines 7a(t'^) of the emitters in an inertial system {x"}, as well as the space- 
time metric in it, are here supposed known 'background' data for any user of the 
positioning system (the world-lines can be pre-determined initially or broadcast in real 
time). Thus, any user who receives only the emission data {t^} is able to compute the 



quantities 74, y*, x appearing in (40). He has to follow the four steps: 

Step 1. Compute the four position vectors of the emitters, 7a = 07a(t^), for the 
received values {r^}. 

Step 2. Choose a reference emitter, say 74, and compute the position vector of the 
three referred emitters = 7a(T") — l^ij'^)- 

Step 3. Compute the configuration scalars VLa = ^e^, the configuration vector 
X = *(ei A 62 A 63), and the configuration bivectors E"- = *{ea+i A ea+2) and H = QaE"". 

Step 4- Choose a transversal vector ^, ' X 7^ 0, and compute y^: = ■^i{^)H. 

At this level, the user has computed the quantities 74, y*, x- But he is also able to 
compute in what of the coordinate regions, Ct, Ci or of the positioning system he is; 
for this, one additional step is sufficient: 

Step 5. Determine the sign of ex^- 
Then, according to their definitions, the user is in Ct, Ci or Cg if this sign of ex^ is > 0, 
= or < respectively. 

Now, suppose that the user is in Q or Cg- Another step allows him to compute the 
orientation e: 
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Step 6. Determine the sign oiu-x for a future-pointing time- like vector u, arbitrarily 
chosen. 

Then, according to proposition [9| e = sgn{u ■ x)- 

We shall call central region of a positioning system the region = Cs U C^. A part 
of theorem [2] may be then stated as follows. 

Proposition 10 Let'-y^ be the position vectors of the known world-line equations 7a(''"'^) 
of the four emitters with respect to an inertial coordinate system. Then the users of the 
central region C'-' of the positioning system can obtain their position in the inertial 
system by computation from their sole emission coordinates {t^}- 

What about the users in Ct, i.e. out of the central region We have seen that 
everywhere e = sgn {y ■ x)- Nevertheless, the null vector y involves not only the initial 
data 7a(t'^) and the received data t"^, but is the solution we are looking for, so that 
the quantity y ■ x cannot be computed before itself. For this reason proposition [9] is 
exclusive: only the users of the central region C*" are able to compute the orientation e 
of the positioning system. 

Can the users in Ct know the orientation of the positioning system? It is possible 
to show that all the users of the coordinate region C, be them in the central region C*^ 
or not, can determine the orientation of the positioning system if, in addition to the 
reception of the {t"^}, they are able to observe the emitters in their celestial sphere. 
But this fact will be analyzed in a forthcoming paper. 



9. Discussion and work in progress 



In this paper we have obtained the coordinate transformation (40) between emission 
coordinates and inertial coordinates in Minkowski space-time. We use the intrinsic 
vector formalism to express the position vector, x = (x"), of every event in the emitter 
coordinate region as a function / of the emitter world-lines 7a(''"'^) : x = /(07a(t'^)) 

This general and compact expression of the coordinate transformation will be a 
powerful tool for subsequent applications. For example, we can particularize it for 
different choices of the emitter world-lines which model specific physical situations. In 
doing so, a previous basic task appears to be convenient in many cases: to write our 
covariant expressions in a 3+1 formalism with respect to an arbitrary inertial observer, 
which is the goal of another work [15]. Moreover, from the expression of the coordinate 
transformation we can easily obtain the components of the metric tensor in emission 
coordinates and we can study the region where the emission coordinates are more 
efficient than the inertial ones [16]. 

In section [6] we have studied the emission conditions in order to distinguish between 
emission and reception or mixed configurations. These results are a necessary tool 
to carry out more in-depth analysis of the domains and co-domains of the emission 
coordinates. This task will be tackled in another paper where we will also study 
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the geometry of the emitter configurations attending to their different causal character, 
and we will analyze how this geometry infiuences the solutions to the null propagation 
equations. Some preliminary results on this question have been presented in [7]. 

The orientation e of a positioning system with respect to an event (see section 
[7]) is a concept which has allowed us to give an explicit expression of the coordinate 
transformation in theorem|2j In section|8]we have pointed out that this orientation e may 
be computed from the emission data {t^} in the so-called central region. Nevertheless, 
out of this region, the determination of e demands an observational method which will 
be analyzed elsewhere. 

The development of all this theoretical work paves the way to the study of more 
realistic, gravitationally infiuenced, positioning systems. For example, those defined 
by emitter world-lines modeling a satellite constellation around the Earth in a weak 
gravitational field. It is worth remarking that for the use and smooth running of 
a positioning system one needs, not only the coordinate transformation to emission 
coordinates but also a good understanding of the domains in the space-time and co- 
domains in the grid. The analysis of the degenerate configuration data, x = 0? ^ind 
of the hypersurfaces of the emission region TZ where the Jacobian je{x) vanishes must 
also be well understood. We already know that the events of vanishing Jacobian are, 
and only are, those for which any user in them can see the four emitters on a circle in 
his celestial sphere [5]. This includes the possibility for the user seeing less than four 
satellites, when some of them are in the shadows of the others. 

Obviously, in realistic situations, the inertial coordinate system considered here 
must give rise to more useful ones. Namely: the International Celestial Reference 
System (ICRS) for the positioning system based on four millisecond pulsars, valid 
for the Solar System, as proposed in [ISj; the Barycentric Celestial Reference System 
(BCRS), to compare the planet trajectories (Earth at least and already) obtained in this 
millisecond pulsar system with the ones obtained by standard astronomical observations; 
the Geocentric Celestial Reference System (GCRS) or related World Geodetic System 
84 (WGS84) or International Earth Reference System (ITRS) for the applications of 
relativistic positioning systems to the Global Navigation Satellite Systems (GNSS). 
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Appendix A. Proof of proposition [l] 



From lemma [2| to prove proposition [T] we must obtain a particular solution of the main 
linear system (10). More precisely, we must obtain the particular solution y^, which is 
orthogonal to a chosen transversal vector ^, ^ ■ X 7^ 0- 
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In order to obtain such a solution in a covariant way we begin by studying a 
regular linear system 0^ ■ 2; = Q^. In this case, the vectors {(Pa} define a base, that 
is, 01 A 02 A 03 A 04 7^ 0. Then, the sole solution to the system takes the expression 
z = Qa (p^, {(p^} being the dual base. More explicitly, we have: 

Lemma 5 The solution z to the regular linear system 

<Pa-z = VLa, a = 1,2,3,4, (A.l) 



is given by 



where 



nA<P^ ■ (A.2) 

= ^e^'''^'' * (0p A 0Q A 0k) , D = *(0iA02A03A04) . (A.3) 



Now, we regularize the main linear system ( 10 ) by adding a new equation as follows. 
Given a vector ^ such that ,^ ■ x 7^ 0, let us consider a linear system of the form (A.l) 
for the unknown z, with 

0A = ea, VLA = VLa, A = 1,2,3; (A.4) 
04 = ^ , ^4 = . (A.5) 

This linear system is regular since: 

D = *(0iA02A03A04) = *(ei AeaAesAO , , 

(A. 6) 

= - * A ei A 62 A 63) = * (ei A 62 A 63) = ■ X 7^ . 

Consequently, the sole solution to this system may be obtained as stated in lemma [5] 
Now f24 = 0, and then we only need to calculate the vectors 0^, a = 1, 2, 3, of the dual 
basis given in (A.3). They take, in this case, the expression: 



3^e»^Q« * (0P A 0Q A (Pr) = 2^e"''^4 * (e, A e, A 
2^e»^^ * (e A A ej = -^e'^''%0 * (e^ A ej = ^tiOE'^ 



(A.7) 



where 



E'^ = *(e,+i A e,+2) • (A.^ 



Then, the solution (A.2) becomes now 

z = -^ tiOH , H = n^E'' . (A.9) 
^ ■ X 

Finally, note that this vector 2; is a solution to the main linear system and it is orthogonal 
to ^. Thus it is the particular solution ?/* that we are looking for. 
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